WIS-99/14 Mar-DPP 



Extraction of nucleon momentum distributions from inclusive 

electron scattering on nuclei 

A.S. Rinat and M.F. Taragin 
Department of Particle Physics, Weizmann Institute of Science, Rehovot 76100, Israel 



(February 9, 2008) 



u . 
< 

Abstract 



> 

\ We address the problem of extracting single-nucleon momentum distribu- 

<0 tions n(jp) from inclusive electron scattering data. A model for these relates 

0\ 

nuclear and nucleon structure functions (SF) through an intermediate SF f PN 
for a nucleus of point-particles. In addition to the asymptotic limit (AL) which 
depends on n(p), f PN contains generally (/-dependent Final State Interactions 
(FSI) parts. In the inverse problem one wishes to separate (/-dependent FSI 

•rH . 

from the AL. In general it suffices to know the structure of the theory, but 
not numerical results. It appears, that in the (/-range of the analyzed electron 
scattering data, FSI are only weakly (/-dependent, making it virtually impos- 
sible to obtain parameters in a free fit of the parametrized components of /. 
Imposing a restriction, we obtain n(p) for Fe and 4 He. 
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The simplest bulk properties of a many-body system in its ground state are the number 
(charge) and momentum distributions of the constituents, or of their centers if the latter 
are composite. Their obvious importance invites measurements, allowing the extraction of 
these observables in a model-independent fashion. This appears feasible for nuclear charge 
distributions, using precision data on elastic electron scattering and on muonic atoms@. The 
situation is different for single-nucleon momentum distribution (MD) n(p). Contrary to the 
case of charge densities, attempts to determine n(p) have met with difficulties, which are 
frequently circumvented by calling in theoretical MrJUl. 

We illustrate the above for exclusive A(e,e'p)(A — l) n and semi-inclusive A(e,e'p)X re- 
actions. An analysis usually starts with the Plane Wave Impulse Approximation (PWIA) 
which expresses yields in terms of single-hole spectral functions. For the two types of re- 
actions these functions are in turn respectively linked to the occupation probability n a of 
a single-particle orbit a and the total momentum distribution n(p) i. However, the PWIA 
result is distorted by Final State interaction (FSI) which as a rule cannot be removed exper- 
imentally. An outstanding exception is the case of near-elastic, high-g electro-dissociation 
of the D, which produced n D (p) i. 

A similar unsatisfactory situation holds for totally inclusive processes A(e, e')X, which 
may also be treated in the PWIA and also require the calculation of FSlSi. Similar problems 
beset indirect attempts to obtain MD. These are based on extracted scaling functions £(q, y) 
which depend on the momentum transfer q and some scaling variable y. In particular 
Gersch- Rodriguez-Smith (GRS)-type of theories^ directly relate lim^oo £(g, y < 0) to the 
MD. Plots of £(q,y < 0) as function of y display coarse-grained scaling, i.e. clustering of 
data points for different q as function of y. Perfect scaling is interpreted as the absence 
of FSI, in principle enabling the extraction of n(p), while imperfect scaling manifests FSI 
0^0. The latter is again difficult to remove experimentally. 

In an alternatively approach one searches for a plateau in £(q, (y) < 0) for binned (y) as 
function of q or of the squared 4-momentum Q = q 2 — v 2 and which is associated with the 
asymptotic limit (AL) for that y, and thus with n(y) BSE!. Both approaches assume that 



FSI terms are recognizable by a clear Q 2 -signature. 

We now proceed and present a different method which proved eminently succesful in the 
extraction of MD from data on the SF (or response) <p(q, y) of liquid 4 He 0^0 and Ne& 
For systems with smooth inter-particle interactions V, SF may be expanded in a series in 

/ A/fx n 

<f>(q,y) = E - F n(v) = <P as (y) + <P odd (q,y) + ^ even (q,y) (ia) 

n>0 ^ ? ' 

F (y) = lim <j>{q,y) = / dppn(p), (lb) 

<P odd (q,y) = E (-) Fn(y) = u^(q)yJ2^ ) y 2n eM-[A {0) y] 2 ) (1c) 



n, odd 



r en (q,y)= E (-) ^d/) = f/ {e) (9)E^ e) ?/ 2nex p(-^ (e) i/] 2 ) ( ld ) 

n,even ^- Q * n 

The above decomposition is in terms of coefficient functions F n (y), even and odd in y, which 
are preceded by factors with?, respectively? even and odd powers of 1/q. The even AL, Eq. 
(|Ta|), depends in a simple fashion on the MD. 

Eqs. (1) instruct how to calculate (f>(q,y) from F n (q,y), which requires as input n(p), V 
and special density matrices^. In practice this is feasible for the AL and for the dominant 
FSI coefficients Fi(y), F 2 (y^^. Alternatively one may wish to use SF data, in order to 
conversely extract n(p). Eq. (|ia]) then serves to express information on the underlying 
dynamics, which one incorporates in parametrizations as in (|lc|), (0)- Algorithms show 
how to separate g-independent AL from, in principle g-dependent, FSI without knowledge 
of computed F n . From fits, including those for the AL, one generates a MD, which appears 
to be in good agreement with accurate calculationsEHl!. 

Two dynamical extensions complicate the above NR approach. The first one arises if 
the bare inter-constituent interaction is not smooth or even singular, which requires the 
replacement of V by V e ff = t q , the t-matrix associated with V. Since V e ff depends on q, 
the GRS g-signature of FSI as in ([lc|), ( p^ will be modified, although for quantum gases 
not drastically^. 

A second complication occurs for molecular systems with additional degrees of freedom, 



which may be excited along with the motion of the molecular centers of mass. Prime 
examples are rotations and vibrations in a gas of H 2 molecules. Again, because one knows 
the required additional dynamics for NR systems, the above complications can be handled, 
both in calculations of the response and in methods to extract n(p). As long as FSI parts 
possess a distinct g-dependence, the algorithms remain applicable!!!. 

We now turn to inclusive electron scattering from nuclei, for which the cross section per 
nucleon reads 



d 2 a eA {E;6,v)/A 2_ /E1 _ n _ A \xM 2 _ a . , + i, Qlt *„ A( ^ 



dTtdv A 



cr M (E; 6, v) 



-Ff(x,Q A ) +t8,n'(9/2)Ff(x,Q A 



(2) 



Q 2 

The inclusive, as well as the Mott cross section <tm, are usually measured as functions 
of E,8,u which are the beam energy, scattering angle and energy loss. The above SF 
Ff 2 (x, Q 2 ) describe the scattering of unpolarized electrons from non-oriented targets and 
are conventionally expressed as functions of the squared 4-momentum Q 2 = q 2 — v 2 and the 
Bjorken variable x = Q 2 jIMv with range < x < A. The strong variation in both cross 
sections in (0) can be exploited by using the tempered, dimensionless ratio 



h A {E- v,B) = —( eA ) ' 1 >' = -pvrFfix, Q 2 ) + tan 2 (^/2)F{ 4 (x, Q 2 ) 



(3) 



2V c7 M (E;u,9) ) Q 2 

In spite of technical complications, due to a finite number of constituents which obey Fermi 
statistics, a description of inclusive scattering of relatively low-energy leptons is not basically 
different from the above NR case. In particular one may still venture to use the notion of 
a potential for the description of iViV collisions in FSI at medium q lab m q < 0.5 GeV. 
However, data for the above kinematics yield only information on MD for restricted p. In 
order to extend that range, one needs considerably larger momentum and energy transfers 
q, v. Those are provided by multi-GeV beam energies, which may excite sub-nucleonic 
degrees of freedom. 

Contrary to NR systems, there is no way to accurately compute SF for nuclei with 
composite constituents. Whatever the approach, it seems evident that with nucleons as 
major constituents, one has to relate nuclear SF in (0) to those of the nucleon Fj? , which 



in practice is the p, n- weighted Fj!' n close to FjP/2. We shall use below such a previously 
formulated relation, which refers only to nucleonic and sub-nucleonic degrees of freedom!!. 
It specifically disregards virtual cloud pions (see ref. [24] []. 

F k A (x,Q>) = J x A ^- k f PN (^Q 2 )F k N ^Q 2 ) (4) 

The above contains SF for a nucleon F{? and f PN for a nucleus of point-particles, which 
has to be computed. Its asymptotic limit depends on the MD of point-nucleons (or of 
the centers of composite ones) and on FSI parts, which account for the distribution of the 
energy-momentum transfer to several core nucleons through iViV collisions. In the kinematic 
range of interest those occur at relativistic momenta q, for which the notion of a local, 
energy-independent potential breaks down. Again one needs V e ff = t(q), which in general 
is off-shell, yet can be parametrized in terms of observable NN scattering data, and which 
permits a calculation of the analog of the SF (p(q, As a first step one replaces there 

the NR GRS-West scaling variable by a relativistic oneil 

Mv ( q 2 \ Mv ( Q 2 \ M . 

y = ~r V ~ m, ) -» VG * -q- v - 2M~v ) = \ in 7 m (1 " x) (5) 



q \ IMvJ q \ IMP J ^l + AM 2 x 2 /Q 2 ) 

Next one establishes the following transition from the SF 0(g, yo) to a relativistic analogc 



r N (x,Q 2 ) = MD{x,Q 2 )<P[q{x J Q 2 ) } y G {x,Q% (6) 

where D is a kinematic factor, which guarantees proper normalization of f PN . The 
theory@§ based on (H) accounts well for the data§§ 



x Eq. @ is easily interpreted in terms of momentum fractions (MF): The MF of a quark in a 
nucleus is the product of MF's of a quark in a nucleon and of a nucleon in a nucleus. SF are the 
probabilities of their occurrence. The Q 2 — > oo limit of a MF equals the observable Bjorken variable 
x = Q 2 /2Mv. Eq. conjectures a similar relation for large finite Q 2 . Similar expressions have 
been derived in perturbation theoryi! with the nucleon off its mass-shell. Fj? in (f4|) relates to free 
nucleons. 
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We now consider the inverse problem, i.e. the extraction of MD from inclusive scattering 
data. We choose the recent Fe data for E = 4.05 GeVi and older 4 He data for E < 
3.6 GevS and determine the reduced cross sections fiA, Eq. (|3|). Its right hand side is 
expressed by means of the nuclear SF (^), which the nucleon SF Fj™. The latter can be 
decomposed as F N = F ,NI (x, Q 2 ) + 5(1 — x)F N, (Q 2 ), containing the nucleon inelastic (NI) 
and nucleon elastic (NE) parts. The former have been parametrizedliS, whereas the latter 
are combinations of known static form factors. 

Eq. (P) relates f PN to 0(g, y G ), with parametrizations ([laD as in the NR case. In 
particular its AL ([TED is generated by a MD, for which we choose a sum of two centered 
gaussians0 

n(p, lk) = (exp[-(p/p!) 2 ] + e exp[-(p/p 2 ) 2 ]) (7) 

In initial runs we found it impossible to obtain a free fit for the FSI part and decided to 
generate starting values for parameters, by fitting the theoretical FSI to (|ic|), (ITJ3). In the 



relevant g-range 0* h (g, 2/g) i n ( P^D appears to have the following striking features: 

i) cj) odd is negligibly small, while the remnant even FSI part A<p even may reach 35% of 
F (y G ) for \y G \ w 0.2 GeV, and up to 15% for \y G \ w 0.4 GeV. 

ii) cj) odd and A<fr even are only weakly g-dependent. 

Both features are properties of an effective interaction of a diffractive nature, for instance 
generated by a nearly imaginary V e ff = t q which interchanges the roles of odd and even parts 
in y G . In the relevant range, its strength Imt g oc a tot happens to be hardly g-dependent, 
which spells difficulties for the separation of the AL from FSI with that property [] Additional 
evidence may be implicit in an analysis by Ciofi et a/H who found that extracted g-dependent 



x Not only NN interactions produce a diffractive a q l '. For example, the atom-atom interaction 
in liquid 4 He has a strong, short-range repulsion!!. This causes damped oscillation in 0"*°* for 
increasing q. However for all, but the largest relevant q, there still is a discernible g-dependenceH. 
The latter property and the relative weakness of FSI, enable the extraction of the MD from the 



scaling functions can be fitted by a g-independent parametrization. This surely is compatible 
with weakly g-dependent FSI. 

We return to attempts to extract n(p) from Fe data, which we divided in 10 (q) bins. 
The lowest two, which correspond to 9 = 15, 23°, resist any fit. This observation is in 
line with the validity of the relation (|j), which predicts deterioration with decreasing Q 2 . 
For increasing Q 2 , data in a given (q) bin often produce excellent fits for hp e with relative 
deviations from data, rarely exceeding 5%, and usually staying under 2%! Nevertheless, the 
corresponding parameters occasionally appear far from their starting values. 

The failure of a free fit forces the imposition of a restriction, e.g. a given central value 
n(0) in the neighbourhood of the theoretical one. Such a choice leads to a well-determined, 
reasonable value for pi, the width of the dominant gaussian in the MD ([F]). In contra 
distinction, the width of the second gaussian appears strongly correlated to its relative 
strength e and the data only approximately determine e 1//3 p2- All fits produce a smooth 
n(p) which are primarily differentiated by their extention, i.e. by the rms momentum. For 
He we could not reproduce the small starting value e = 0.003 and just put it to 0. Table 
I assembles starting values and fitted parameters. The shaded area in Fig. 1 shows the 
extremes for n Fe (p), corresponding to the entries in the Table. The drawn line represents 
the MD for the given starting values. Fig. 2 for He presents only one fit. 

We summarize. Encouraged by the succesful extraction of single-atom momentum distri- 
butions from cross sections for inclusive scattering of neutrons from mono-atomic quantum 
gases, we attempted such an extraction fron nuclei. The feasibility in the case of NR sys- 
tems rests on knowledge of the atom-atom interaction in the above cases, and of additional 
rotation-vibration dynamics in the case of di-atomic molecules, etc. This knowledge en- 



dominating AL. The most extreme example is that of hard-core interactions, which produce con- 
stant a* *. At least part of the FSI are strictly g-independentB The above implies that in the 
large Q 2 -plateau of a plot of £,(q,yG) nearly g-independent FSI 'contaminate' the AL there. 
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ables parametrizations of the structure function, incorporating the above information. Fits 
to data, enable a separation of g-dependent FSI parts from the g-independent asymptotic 
limit. That AL has a simple dependence on the desired momentum distribution, which is 
extracted in a truly model-independent manner. 

The nuclear case is in many respects much more complex. For one, the composite nature 
of nucleons cannot be accurately ?order? described as in the molecular case. Consequently, 
one has to invoke some model dependence, e.g. the one implicit in ([§). Moreover, there 
is a real stumbling block in the shape of FSI, which are barely g-dependent. Consequently 
a free fit to the data does not lead to fits, smooth in the parameters. Only restrained fits 
determine characteristics on the momentum distributions. 

The authors acknowledge stimulating discussions with many colleagues on the subject 
matter. Special thanks go to Tim Shoppa and in particular to Byron Jennings, who partic- 
ipated in initial stages of this research. 
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Figure captions 

Fig. 1 Momentum distribution n Fe (p) from fits of d 2 a III for parameters in Table I. 
Shaded areas correspond to entries in Table I with rrii corresponding to different fixed n(0). 
The drawn curve is for the starting values of the parameters^. 

Fig. 2 Momentum distribution n He (p) from fits of d 2 a for parameters in Table I. The 
drawn curve is for the starting values of the parameters^. 

Table I 
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Fe 


He 


Parameters 


Starting values 


extracted 


extracted 


n(0)(in fm -3 ) 


69.8 
76.8 
58.3 


80.51 


Pi (in fm^ 1 ) 


0.802 


0.837± 0.011 
0.813± 0.009 
0.797± 0.041 


0.821 


e 1 / 3 p 2 (in fm -1 ) 


0.432 


0.406±0.060 
0.365±0.060 
0.268±0.036 




p 2 (in fm -1 ) 


1.390 


1.177±0.136 
1.371±0.206 
1.235±0.211 




e 


0.03 


0.041±0.020 
0.022±0.009 
0.020±0.010 


0.0 


\/(P 2 )n(p)(in fm" 1 ) 


1.105 


1.108 ±0.035 
1.105±0.037 
1.013±0.021 


1.005 



Extracted parameters for n(p) of the form Eq. (|6]). For three fixed values of n Fe (0) 
are entered the width of the dominant Gaussian pi, the combination e 1//3 p2 of the relative 
strength e and width of the second one, as well as each separately, and the root mean square 
momentum. For He we present only the fit for e=0. 
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